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Following our previous work, a complete classical solution of the CGHS model in Hamiltonian 
formulation in new variables is given. We preform a series of analyses and transformations to 
get to the CGHS Hamiltonian in new variables from a generic class of two dimensional dilatonic 
gravitational systems coupled to matter. This gives us a second class system, a total Hamiltonian 
consisting of a Hamiltonian constraint, a diffeomorphism constraint and two second class constraints. 
We calculate the Dirac brackets, bring them to a standard form similar to the Poisson brackets 
by introducing a new variable. Then by rescaling lapse and shift, the Hamiltonian constraint is 
transformed into a form where it has an strong Abelian algebra with itself. This property holds 
both in vacuum case and in case with matter coupling. Then for each of the vacuum and the coupled- 

*vj ' to-matter cases, we preform two gauge fixings, one set for each case, and solve the classical system 

completely in both cases. The gauge fixing of the case coupled to matter is done by implementing 
a method based on canonical transformation to a new set of variables and leads to a local true 
Hamiltonian. 

We also show that our formalism is consistent with the original CGHS paper by showing that the 
equations of motion are the same in both cases. Finally we derive the relevant surface term of the 

C*J • model. 
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I. INTRODUCTION 

The two dimensional gravitational systems with black hole solutions, specially the CGHS model [l|, have proven 
to be a very good test bench to try out various ideas about quantum gravity. There is an extensive study of these 
systems in the literature (for a short list of references see [2|t8| and the references within them) and there have 
been numerous attempts to understand some quantum phenomena, such as black hole evaporation, information loss 
and the asymptotic fate of spacetime, using their black hole solutions. There is also the important question of 
whether quantum gravity, specially loop quantum gravity, eliminates the singularity? Loop quantum gravity has been 
progressing on the issue of addressing how singularities are affected by quantum gravity. Examples are replacing 
the Big Bang singularity by a Big Bounce in homogenous models [9] and some form of singularity resolution for the 
Schwarzschild black holes in spherically symmetric midi-superspaces [10]. But with all these attempts, the questions 
surrounding these issues has not been answered in a satisfactory way. 

The purpose of this paper is twofold. On one hand, we address the formulation of the CGHS model in a manner 
that is suitable for applying loop quantum gravity techniques. We start by deriving the Ashtekar-like variables for the 
model and writing its Hamiltonian in terms of those variables. One of the differences of our method with many other 
works that can be found in the literature is that this formulation is pursued without a conformal transformation. 
This is important in the sense that one is working with variables that have direct geometrical meaning so there is 
no need to turn everything back to their original directly-geometric form at the end. It also makes it easier to read 
the physical implications off of the theory. Another advantage of our formulation is that it gives us a Hamiltonian 
constraint that commutes with itself and hence the algebra of the constraints becomes a Lie algebra, since now we 
have structure constants instead of structure functions. Thus this formulation of the system appears to be suitable 
for applying loop quantum gravity techniques. 

The second purpose of the paper is to solve the classical system coupled to matter completely classically by some 
choice of gauge fixing such that it results in a local true Hamiltonian for the system. Local here means that the lapse 
will not be an integral of canonical variables and hence the Hamiltonian will not be an integral of an integral. This 
way the equations of motion will also be local. The locality is important because, among other things, it is obviously 
much harder to try to quantize a non-local theory. There have been some previous studies [1 ll — tl3l | which did not 
lead to a Hamiltonian that was the spatial integral of a local density, thus leading to non-local equations of motion. 
This in turn leads to difficulties upon trying to quantize the system. However, there have been some recent works 
|14lfla| leading to local true Hamiltonian for 3+1 spherically symmetric case which we follow in order to derive such 
a Hamiltonian for the CGHS model in new variables. 

The structure of this paper is as follows: in section [Hi we rewrite a generic action for two dimensional gravitational 
systems coupled to a dilaton field and a scalar matter field (which includes the CGHS model as a special case) in a 
suitable way to include cases with or without kinetic term for the dilaton field in one action. In section Hill the previous 
Lagrangian will be written in tetrad variables and is transformed into a Hamiltonian by a Legendre transformation. 
In section HVl we focus just on the CGHS model, derive the new Ashtekar-like variables for it similar to the 3+1 
case and write its Hamiltonian in those variables. Sections [V] and IVII are dedicated to the implementation of the 
Dirac procedure for a second class system in our Hamiltonian formalism. In section fVIIi the Hamiltonian constraint 
is written in a special form such that it has an strong Abelian algebra with itself. Section IVIIII is where equations 
of motions are derived for the Hamiltonian system. In sections IIXI and [Xj both vacuum and coupled-to- matter cases 
are gauge fixed and a local true Hamiltonian is derived for the case coupled to matter, section IXII is dedicated to 
show that our Hamiltonian formalism is equivalent to the original Lagrangian formalism by comparing the equations 
of motion in both cases. Finally, in section IXII1 the boundary term for the formalism is derived and is compared to 
the standard boundary term of the CGHS model. 

II. GENERIC ACTION IN METRIC VARIABLES 

It has been shown @, Hs, LLZ] that the most general diffeomorphism invariant action yielding second order differential 
equations for the metric g and a scalar (dilaton) field <& in two dimensions i£j 



Vdii 



J d 2 x^\g~\ (D($)R(g) + l -g ab d a $>d b <5> + V{$)) . (1) 



1 The authors in Q argue that the most general form is in fact S = f d 2 x^/—\g\ I D(&)R(g) + V ( (V$) 2 , <£ ) ) . 



In fact two of the most important models of gravitational systems are examples of this action. The first case, the 4D 
spherically symmetric model, is just the Einstein-Hilbert action. Minimally coupled to matter it is 

Ssph = ^ f dWHlfa - i- f d 4 xy/^W\g ab d a fd b f, (2) 

which by using spherically symmetric ansatz 

ds 2 =g liV dx^dx u + $ 2 (d6> 2 + sin 2 ^)^ 2 ), fi, v = 0, 1 (3) 

and integrating over and (/>, becomes 



Ssphcr - J d 2 x^\g~\ Q$ 2 i?( 5 ) + \g ah d a <S>d b <5> + \)-lf d 2 x^\g-\$ 2 g ab d a fd b f. 



(4) 



In (|2]), x°, x , and 4> are some coordinates adapted to the spherical symmetry and g^ v is the metric on the x , x 1 
plane. If Jp- ^ 0, it is possible to rescale $ to be identical to x 1 . It is easily seen that (the gravitational part of) (0} 
is an example of ([T]) and for this case, $ mimics the dilation field. 

The second case is the CGHS model [1] whose action minimally coupled to matter is 

Scghs = J d 2 x^\g~\e- 2 v (R + ±g ab d aV d W + 4A 2 ) - i j d 2 x^\g~\g ab d a fd b f, (5) 

in which tp (not to be confused with the coordinate (f> in spherically symmetric model) corresponds to the dilaton field 
and A 2 is the cosmological constant. By a redefinition of the dilation field 

$ = 2%/2e^ v , (6) 

one gets 

Scghs = J d 2 x^\g~\ |l$ 2 i?+ I^ b a a $a 6 $ + i$ 2 A 2 | -\fd 2 x^g~\g ab d a fd b f (7) 

whose gravitation part can be seen to be an example of ([1]). So from these two cases we can infer that the minimally 
coupled matter part corresponding to ([lj, can be written as 



S m = - j d 2 x^\g~\W(<S>)g ab d a fd b f. (8) 

Thus the full general action with minimal coupling to matter will become 

S = J d 2 x^\g~\ (D($)R(g) + \g ab d a <S>d b <5> + V{$)\ - j d 2 x^\g~\W(<P)g ab d a fd b f. (9) 

In some cases, it is desirable to eliminate the kinetic term g ab d a Qd b &. This can be achieved by a conformal transfor- 
mation 

g ab = n 2 ($)g ab , (10) 

with 

and C being a constant of integration. In this case the general action becomes 

S = J d 2 x^\g~\ (D($)R(g) + Q- 2 ^)V(^)) - J d 2 x^\g~\W^)~g ab d a fd b f. (12) 

In order to keep the generality, we can combine ([9]) and (fl"2"T) into the following form 

Sx+x = j ' d 2 x^\g~\lY(<$>)R+^Zg ab d a $>d b <5> + V($)\ - J d 2 x^\g~\W($)g ab d a fd b f, 



(13) 



(14) 



where we have introduced the variable Z = 0,1, which plays the rule of a "switch" that turns the kinetic term on or 
off. The gravitational and matter Lagrangian densities are now 

L % = ^\g~\ h($)R + \zg ab d a <$>d b <$> + V($)\ , 

L m = -^\g~\W($)g ab d a fd b f. (15) 

It is worth noting that the action (fT3|) is similar to the action of f(R) gravity theories with V($) being the potential 
of the dilaton field. In the CGHS case, this potential is related to the cosmological constant. So the cosmological 
constant in this case may be seen as dynamical and coming from properties of the dilaton field rather than just being 
there. 

III. TETRAD FORMULATION 

In order to get to the new variables for these cases, one needs to first write the theory in tetrad formulation in 
which 

g a b = ViJe I a e J b (16) 

and rju is the Minkowski metric, e 1 a are the tetrads and /, J are the internal indices while a, b are the abstract ones. 
The curvature can be written in terms of the curvature of the spin connection and ultimately in terms of the spin 
connection u! a IJ itself as 



R = Rab e je 



= (2d [a LJ a] IJ + [u a ,Ub\ IJ ) e Q /e" j, (17) 

where [, ] stands for the Lie commutator in the Lorentz Lie algebra and the indices /, J take value in this algebra. 
Since the spin connection is antisymmetric in /, J, we can write it as 

uj a IJ = to a e u . (18) 

This way, the curvature (fTT|) becomes 

R = (2d [a Lo a] e IJ + ^ [aW6] 77 /J ) e a je b j 

= 2d [a uj a] e IJ e a ie b j, (19) 

where we have used the following fact about the Lie commutator in this case 

[u ai uj b ] IJ = uj a e I K u! b e KJ - OJbe 1 K ui a e KJ 

= u [a u b] e 1 ^' 1 (20) 

and also the fact that r] IJ e a ie b j is symmetric in a, b while dJ[ a u> b ] is antisymmetric and thus 



Wr, 



[a uj b] T 1 JJ e a I e b j = Q. (21) 

We also would like to add the torsion free condition (contracted by a Lagrange multiplier) to the action ([LI)) . This 
condition reads 

= de 1 + e 1 jlu A e J = 2d [a e b] I + 2e T juj [a e b] j . (22) 

Since this is a 2-form, we need to contract it with e ab to get an scalar to be able to add it to the Lagrangian density. 
Doing so, contracting it with a Lagrange multiplier — Xi and substituting everything in (|14[) yields 

L g = -2X / e a6 (9 [a e h] / + e 1 jLu [a e b] J ) + 2Y8 [a uj b] ee IJ e a ie b j + ^Zr) IJ ee I a e J b d a $d b $ + eV. (23) 

We can write e ab in terms of e IJ , eJ and its determinant e = det(e a 7 ) as 

e ab = -ee/ejV' 7 . (24) 



If we integrate by parts in the first term in Lagrangian density (f2"3"f to bring the partial derivative to act on Xi, and 
then use the above result for e ab in both first two terms of this Lagrangian, the pure gravitational Lagrangian density 
can be written as 



L g =e (-2d a {Xj)e K a e KI - 2Xje Ia uj a + 2Yd a oj b e IJ e I a e J b + ^ V J e/ Q e/d a $c> 6 $ + v) . 



From (|f 5| . the matter Lagrangian density can also simply be written as 

L m = -W V IJ e ei a ej b d a fd b f. 



(25) 



(26) 



Next, we decompose the Lagrangian by ADM method and perform a Legendre transformation to get to the Hamil- 
tonian. Most of the details needed for these steps have been already discussed in our previous work [l8| . So we just 
mention the result here. The generic Hamiltonian then will become. 



H =N 



2 w\ { * xl) 



Z $'2 _ J|_ + 



Pi 

{ 

2Wf 2 



2 w\ (* x2 y 



2 wr* xl ~ 2 wr* x2 



n 



\p 



N 1 



Z\P\ 



2W\P\ 



^ 



P$$' + P f f + Pi (*X 1 )' + P 2 (*X 2 )' - P^X 2 - P 2 uji*X 1 



+ u Ja [P 1 *X 2 + P 2 *X 1 -(2Y)']. 
Here TV is lapse, TV 1 is the (one dimensional) shift vector, the momenta are 



Pr 



P« 



OL 



d*X I 
dL 

dL _ Z^/q 



2y/qru, 



= 2Y, 



9$ 
dL 



p f =— ' = - 
df 



N 
2Wyfq 



$ 



N 



(n 1 ®' 



(27) 

(28) 
(29) 
(30) 
(31) 



n/'s are normals to the spatial hypersurfaces, *X ! = e IJ Xj, q is the determinant of the induced metric on the spatial 
hypersurfaces, \P\ = \J —i~i IJ PiPj is the norm of Pf (which is a timelike vector, hence the negative sign) and the 
prime represents partial derivative with respect to the spatial coordinate x 1 . From the definitions of momenta (f2"9")l 
and (pJU)) . one can see that if Z — 1, i.e. cases with kinetic term present (without conformal transformation), then <I> is 
a canonical variable and Pj> is its momentum q In this case, since Y involves $, equation (|29[) will be a new primary 
constraint. This would not happen in the cases without kinetic term (with conformal transformation) like standard 
spherically symmetric case in Ashtekar variables, since there, <I> is not a canonical variable although equation (|2"9"|) is 
still valid. Thus in the CGHS case, we should add (|2^|) to the general Hamiltonian (|2"T|) to get the total Hamiltonian. 

IV. THE CGHS HAMILTONIAN IN NEW VARIABLES 

From now on we focus only on the CGHS case and hence we substitute the explicit forms of Y($), V(<&) and W / ($) 
in the previous generic Hamiltonian and also add (|29|) as a new primary constraint to it to get 



H -Ml 2P2 Ft *Y 1 J- 2Pl f) *Y 2 2Pl ,,*Fl 2P2 ,,*Y 2 

H-N\—d x X +—d 1 X -— Wl X -— Ul X 



\P\ 



A 2 $ 2 



2 In fact this can be seen more clearly in the Lagrangian which leads to this Hamiltonian. There, all the terms containing <t> are multiplied 
byZ. 



$' 2 n , (/') 2 , p f 



r\ r r ' p _ 

f) 



N 1 (Pxd^X 1 + P 2 di*X 2 - P 2 *X 1 uj 1 - Pi*XV + &P$, + f'P- 



+ uj (p 1 *X 2 + P 2 *X 1 - Q$ 2 

+ m(p uj -U> 2 \, (32) 

where M is a Lagrange multiplier. In order to transform to the Ashtekar variables and following a similar pattern as for 
Ashtekar variables in the 3+1 model [18j, we introduce the following new momenta with a canonical transformation: 

P u =E X , (33) 

\P\ =2E V , (34) 

P 1 =2cosh(r 1 )E ip 1 (35) 

P 2 =2smh{r))E v , (36) 

where E x , E v and r\ are the new momenta. This gives us the generating function 

F(q, P) = 2*X 1 cosh(?7)P v + 2*X 2 sinh^E^ + uj 1 E x + $P<t. + fPf- (37) 

Using F(q,P), we can find the new canonical variables as 

BF 

Q v =— = 2*A 1 sinh(77)£ v + 2*A 2 cosh^)^, (38) 

OF 

K v =-— = 2*X 1 ooshfa) + 2*A 2 sinhfa), (39) 

A.~=^, (40) 

where Q n , K v and A x correspond to 77, E v and E x respectively. From the above equations, we can find *X X , *X 2 
and wi as 

^^(^coshW-^M), (41) 

^ 2 =4(^sinh(,)-^M), (42) 

wi =^x- (43) 



In order to write the Hamiltonian density (|32| in these new variables, we substitute (|33 j) - (|36|) and (|4T[) -(|43 |) in the 
total Hamiltonian (|32p . and then make a field redefinition 

A X = K X - 7/ (44) 

to get rid of r\ in the Hamiltonian and get 

\ EV EV 2 2 v 2EV 2E* 2E"P 2E^ I 

+ N 1 (E*K' V - Q V K X + $'P* + f'Pf) 

+ c (q ri Q$ 2 ) j + M (e x i$ 2 ) . (45) 

We can see from here that the total Hamiltonian is just the sum of four constraints as is expected for a totally 
constrained system. The first constraint multiplied by the lapse function N is the Hamiltonian constraint. The one 
that is multiplied by the shift vector A 1 is the diffeomorphism constraint. The constraint that is multiplied by u>q is 



the Gauss constraint and the last one is the one we got from the definition of the momentum P u . Solving the Gauss 
constraint in the above Hamiltonian and substituting the resultant Q v from it back into the Hamiltonian yields 

'**" **'&" 1™, 2 * 2 „ v Pi , (/') 2 , P f\ 



H =N — ttt^t- ~ ^E V \ 2 <S> 2 - K V K X - 



2Ev 2E^ 2 2 v 2Ef 2Ev 2Ev> I 

+ N 1 (e^K^ - -$><f>'K x + $'P$ + f'Pfj + M (e x - -$ 2 j . (46) 

Since we now know our canonical variables and momenta, we can write their Poisson brackets as 

{K x (x),E*(y)}=6(x-y), (47) 

{K v (x) t E*>{ y )}=6(x-v), (48) 

{$(x),P*(y)} = 6(x-y), (49) 

{f(x),P f (y)} = S(x-y), (50) 

and we have not written the Poisson bracket of (Q v ,v) P a i r because they no longer appear in the Hamiltonian. The 
rest of the Poisson brackets are strongly zero. 

V. THE CONSISTENCY CONDITIONS ON CONSTRAINTS 

Following the Dirac procedure, we should check the preservation of the constraints to see if there are any new 
secondary constraints and/or to find the value of the Lagrange multipliers in terms of canonical variables. This means 
that the constraints C, being the constants of motion, should remain weakly vanishing during the evolution 

C = {C,P}«0, (51) 

where « represents weak inequality. The Poisson brackets of the Hamiltonian and diffeomorphism constraints with 
H vanishes weakly. Let's check the consistency of the constraint 

H = E X -U> 2 . (52) 

For this, we need (|47[) and (|49[) . The preservation condition of /j, constraint leads to a new, and by definition secondary, 
constraint which we call a: 

1 Pf$ , s 

A = {M,#}«0^a = ir ¥ , + --| r «0. (53) 

We also need to check the preservation of the new a constraint. This leads to a relation between the Lagrange 
multipliers N, N 1 and M (and canonical variables). Finding M from this relation and substituting it into the total 
Hamiltonian (I46[) yields 

h = n(- K,jc r - 2 "' EV ' EX + ^.A_V A v i ***. ■ 2P * K * EX 



t-ip-'-^x 



$Ev> 2 <&Ev 2Ef 2 " 2 Ef <£>Ev 

2E X P 2 2<S>' 2 E X $' 2 2E x f' 2 ' 
<I>2£V $2£y 2E^p & 2 Ev 

+N 1 f - )-§&K x + *'P* + f'P f + E X K' X + E V K' V - ^K) ■ (54) 

Next step is to check if the constraints are first class or second class. Calculating the Poisson brackets of constraints 
among themselves shows that ijl and a are second class and do not commute with each other. In other words, their 
Poisson bracket with each other does not vanish weakly. Since there are second class constraints in the theory, now 
we should abandon the Poisson bracket and move on to the Dirac bracket and also put the second class constraint 
strongly equal to zero and eliminate some of the variables in term of others. By doing this, we can get rid of the 
(<&, P$) pair in the Hamiltonian as can be seen below. Equating both the fi and a constraints strongly to zero yields 

fjt = => $ = 2Ve?, (55) 

a = 0=*>P$ = -^L. (56) 

Ve7 



Substituting these in the total Hamiltonian (|54l) yields 

H =N (-K V K X „E^1^^L + E1^^1_ WET* + ^ + ±ll) 
I <p EV 2 2EVE X Ev 2 E x 2 E^ 2 Ef J 

+ N 1 \-K x E" + f'Pf - y Ex + ^ g; j ■ (57) 

where now we are only left with a Hamiltonian and a diffcomorphism constraint. 

VI. DIRAC BRACKET AND THE ALGEBRA OF CANONICAL VARIABLES 

In order to switch to the Dirac bracket, we need to find the general form of the Dirac bracket for our theory. For a 
field theory (where the variables have continuous indices), the Dirac bracket is 

{A(x), B(y)} D = {A(x), B(y)} - J dw J dz ({A(x), Xp (w)}C^(w, z){ Xa (z),B(y)}) , (58) 

where the {, }u refers to the Dirac bracket, %'s are the second class constraints and C p,T (w, z) are the elements of the 
inverse of the matrix of the Poisson brackets between the p'th and cr'th second class constraints 

C pa (w, z) = {Xp(w),Xa(z)}. (59) 

In our model, there are only two second class constraints, fi and a. Thus the matrix of the Poisson brackets of the 
second class constraints will be 

r-r (* ,,\- ({vW'Kv)} Mx),a{y)} 
pa[ ' V) - Ua(x), /*(»)} {a(x), a(y)} 



{n(x),a(y)}\ , . 

{a(x), /*(»)} ■ (W) 



To compute the elements of this matrix we use (|5"2")) and ([531 along with (|49p to get 

1 $(y) 



1<%) 2 
The same method of computations gives 



{$(x) 2 ,F*(y)} 

5{x-y). (61) 



{Mx),a( y )} = i|^^-y). (62) 

To calculate the elements of C _1 , we use the property CC _1 = 1, or in terms of their elements 

j C pa (x, z)C^(z, y)dz = 8/8{x - y), (63) 

which yields 

C- 1 = C»°{x,y)=[ AE , (X) ^)S(x-y). (64) 

Using this and (|55|) . the general form of the Dirac bracket for our theory becomes 

{A(x), B(y)} D ={A(x),B(y)} + J dw J dz (^{A(x)^(w)}^^-S(w - z){a(z), B(y)} 

- j dw J dz l{A(x), a(w)}-^-j-5(w - z){u(z), B(y)}) . (65) 



If we use this formula and the Poisson brackets (|4"7]) - ([5T)1) . we can find the Dirac brackets of the canonical variables 
between each other as 

{K x (x),E*(y)} D = {K v (x),E*(y)} D = {f(x),P f (y)} D = 5(x-y), (66) 

{K x (x),K tp (y)} D = ^S(x-y), (67) 

E v 
{K x ,E">} D = -—5(x-y), (68) 

{E X ,K V } D = {E*,E*} D = {f,U} D = {P f ,»} D = 0, (69) 

where ■ means everything except Pf and the • means everything except /. 

Looking at the above Dirac brackets, we can make an important observation: by introducing a new variable 

E^K 
U* = K x + f^, (70) 

the Dirac brackets (f6"6"]) - ([6"9"|) can be brought to the standard from 

{U x (x),E*(y)} D = {K v (x),E*{y)} D = {f(x),P f (y)} D = 5(x-y), (71) 

with other brackets being zero. This is an important step since it makes Dirac brackets look like the standard form 
of the Poisson brackets of canonical pairs. Using (f70|) . the total Hamiltonian (|57f becomes 

H =N {-K V U X _E^^l_E^ E^ 1 J^ _ + ^ + lll) 

\ v * EV 2 2E^E X E* 2 #" IE* 2E*\ 

+ N 1 {-U x E xl + f'P f + E*K' V ) . (72) 

VII. TRANSFORMING THE VACUUM HAMILTONIAN CONSTRAINT INTO A TOTAL 

DERIVATIVE 

We can omit U x in the Hamiltonian constraint in (|72l by a redefinition of the shift 



"* = Nl + ^f • (73) 



Substituting this into (|T2"j) gives 



/ ^^ 1 ^'2 £*« j^ E*K V K' 
\ E* 2 2E*E X E* 2 E x E xl 

-W^-2E*E X *+ X - P 1+ X J- 
E xi ^Ef 2E* , 

+ N 1 (-U x E xl + f'P f + E V K' V ) . (74) 



Now redefining the lapse in the above total Hamiltonian as 



E 
will yield 



N = N—— (75) 



— ( d (1 E x ' 2 onia;x2 IK%\ fPjKf 1 P]E XI 1 E x 'f' 



H T =N — - - ^ - 2£ a: A 



dx \2 S^ 2 ^ 2 E x I E X E* 2 E* 2 E X 2 E* 2 E X I 

N 1 {-U x E xl + f'P f + E*K') . (76) 



The Hamiltonian constraint above 



d I E*» 1K*\ f'P f K v iFfE" lE*>f* 

dx \ 2 E* 2 E X 2 E x j E X E* 2 E* 2 E X 2 E* 2 E X ' V ' 



10 

has a remarkable property: its complete form or its vacuum form (/ = = Pf) both have a strong Abelian algebra 
with itself, namely 



{H(x),H(y)} D =0, (78) 



H(x) 



My) 

/=0,P/=0 



=0. (79) 

f=0,P f =Q) D 



This way the algebra of the constraints becomes a Lie algebra since now we have structure constants instead of 
structure functions. 



VIII. EQUATIONS OF MOTION 

We are going to write the equations of motion in two equivalent cases, before and after rescaling lapse and shift 
and introducing U x . The reason is that to compare the original CGHS equations of motion with our formulation, it 
is much easier to use the Hamiltonian equations of motion before rescaling lapse and shift. This is because in order 
to be able to compare equations of motion between our formulation and the original Lagrangian formulation, we need 
to find the second order equations of motion from the Hamiltonian equations of motion. This is achieved much easier 
using explicit form of the Hamiltonian before rescaling. 

A. Before introducing U x and rescaling TV and TV 1 

Here we can use the Hamiltonian (|57|) and the Dirac brackets (|M| -(|69 p to find the the equations of motion, 
F = {F(x), J dyH(y)} D , for any function of the phase space F. Using these, one gets for the canonical pairs 



K x 


( K x Ku 1 f' 2 1 Pf E x " E X 'E^' 
VI 1 


E x,2 


x ' I e x ' 2E x Ef ' 2E X E* E x Ev> E X E<? 2 


E x2 Ef 




/,rlrs\> EV ' N ' N " 
V x ) E <p2 fly ' 




E x 
K 


=NK tp + N 1 E x/ , 






*' I 2E^ 2 2E 1 ? 2 ' 2~ "' 2E X ' 2E X E^ 2 J 

N'E XI , , 
l l N 1 K 






1 V u>2 ' ^ A V 





(80) 

(81) 



(82) 

E ip =NK x + N 1 'E ip + N 1 E V ', (83) 

MP 

f=^T+N l f, (84) 

P f = ( N ^+ Nlp f) ■ ( 85 ) 



B. After introducing U x and rescaling N and TV 1 

In this case, we can use the total Hamiltonian (|76[) and the brackets (|7ip to calculate the equations of motion, 
F = (F(x), J dyHxiy)} n - For the canonical pairs we get: 

rj J KyPff f'f" P f P f , E-'(P 2 + f 2 ) \ 

I E x2 Ev E x Ef2 E X E^ E X E^ J 

. f J 1 P /+/ /2 2 E X 'E V ' ^K | E*" 1 E x ' 2 \ 

\ 2 E x Ef 2 E x Ef3 2E x2 E X E^ 2 2E x2 E^ 2 ) 
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AT 



E x =N l E x/ , 



E x 



E x Ef 2 



+ {N 1 U X )'. 



K V =N\ 



EXI (g + / ,2 ) /'P/x; 



£v=iV 



/=iV 



Pff 



E x E<f> 3 



E x Ef 2 



E x Ef E x K ' 

P f E x ' K,J' 



Pt 



N 



f J 

E x Ef> 2 E x Ef 

Wf K v P f 



E X E^ 2 E x Ev 



+ N l f, 



■N l P f 



(86) 

(87) 



rpx/2 

■ N' h N 1 !^' 



(89) 
(90) 

(91) 



IX. GAUGE FIXING THE VACUUM THEORY 



Using (|76|) . the vacuum Hamiltonian (with f = Pf = 0) after an integration by parts can be written as 



H =N'H d + N l V 



=N' I E 



2E'-p 2 E 



1 K 2 \ _i 
- 2E X \ 2 - --£ + N 1 {-U X E XI + E*K' V ) , 



(92) 



where the prime sign on N can be taken to mean that we have a new lapse after integration by parts or to show 
that this lapse is the derivative of the previous lapse TV with respect to x. To completely solve the vacuum classical 
theory, one can use two gauge fixings. We choose the first gauge fixing with an explicit coordinate dependence as 



Xi=E*(x) 



-2\x 



0. 



(93) 



/■ 



Since this gauge fixing does not contain any explicit time dependence, its preservation condition in time gives 

Xi(x) wO (94) 

>l!i{Xilx),H Q (y)} D nO (95) 

N\x)E'(x) swO, (96) 

which implies that the shift (not the original shift but the rescaled one) vanishes, 

N 1 = 0. (97) 

The motivation for the gauge fixing (|93|) is the following: the coordinate transformations 

e \(t+x) 



- + 



A ' 

e -\(t-x) 

X : 



(98) 
(99) 



will lead us to the formulation of the theory in the conformal gauge in the original CGHS paper [lj, where (x~,x + ) 
are the null coordinates in that gauge. Now if we use these transformations in (|93|). the form of E x will becomes 
exactly the same as the form it gets in the vacuum case in the conformal gauge in [l|. 
We choose our next gauge fixing as 



X 2 = E*(x) - 1 » 0. 



(100) 



The reason behind this gauge fixing is the observation that in the conformal gauge, E v is the only independent 
metric component and upon using the coordinate transformations (|98|l and (|99p . the vacuum metric becomes flat (i.e. 



E ip =l). Preserving the \i constraint gives 



which yields 
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X2 (x) wO (101) 

d V { X 2{x),H 1 {y)} D ^ (102) 

N'(x)K v (x)e 2Xx wO (103) 



TV (x) = 0. (104) 



This is not an issue since this is the derivative of the transformed original lapse. Up to now, we have N , TV , E x 
and E v explicitly. Now we can solve for K v from the weakly vanishing of the Hamiltonian constraint in (|92p. This, 
together with JS3]) and (|TDD|) yields 

i^ = 0. (105) 

Finally U x can be found from the weakly vanishing the diffeomorphism constraint as 

U * = "V"' (106) 

which upon substituting the relevant values of the variables from above gives 

U x = 0. (107) 

This way, the vacuum case is completely solved classically. 

A. The original lapse and shift 

Using the results we just obtained for the vacuum case, we can express the original lapse and shift in (|72|) in terms 
of the ones we obtained using the gauge fixings. Using (|75)l . (|104jl . ([93| and (|100|1 we get 

N' =0 (108) 

E V E X \' 



which implies 



for which we can simply choose 



Using ([73]), dm, dH and JT05]) we find 



iV' =0, (110) 



N = g(t) + C, (111) 



N=l. (112) 



TV 1 =0 (113) 

A^) + ^=0 (H4) 

TV 1 ^) =0. (115) 
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X. GAUGE FIXING THE CASE COUPLED TO MATTER 

The above gauge fixings which we used for the vacuum case will not work for the case with matter field. Among 
other reasons, one relatively obvious reason is that the two gauge fixings (|93l) and (|100l) . are related to the form that 
E x and E v take in the vacuum case. Also we can not turn the whole Hamiltonian constraint into a total derivative 
and therefore if we get a vanishing lapse, it would be the original lapse itself that vanishes not its derivative. By 
inspecting the total Hamiltonian (|75|) . one can see that perhaps a good choice for the first gauge fixing is 

Ci = E x - h(x) « 0. (116) 

with h(x) an arbitrary function of x coordinate. Since the only nonvanishing Dirac bracket between £i and Ht in 
(|76|l comes from the term in the diffeomorphism constraint containing U x , and since U x appears without derivatives 
there, the preservation of £i will give us an algebraic equation for N . In more precise way we have 

(i = Ui(x),JdyH T (y)X^0 (117) 

dyN 1 (y){E x (x),U x ( y )}E x '(y)^0 (118) 

dyN 1 (y)E x, (y)S(x-y)^0 (119) 

n\x)E x/ (x)^0, (120) 



which means that the shift vanishes, 



^=0. (121) 



Now, from ([76]) . (11161) and Q12ip . the partially gauge fixed total Hamiltonian will become 

(d(lh' 2 2 lK 2 \ f PfKip 1 P]h> 1 h'f' 2 \ 

HTF - N [dx- \2hE^ ~ 2hX ~ 2—) ~h&r + 2hE^ + 2hE^j ' {U2) 

For the second gauge fixing, we follow the procedure suggested in [1J, [l5| . The basic idea is to make a canonical 
transformation and define a true Hamiltonian that gives the correct equations of motion. To start, we identify the 
terms inside the total derivative of the Hamiltonian constraint in (|122|) as our new canonical variable X, 

1 h' 2 1 K 2 

X = L-^—-2h\ 2 ~^. (123) 

2hE<fi 2 2 h y ' 

This suggests that we can find a generating function of the third kind as 

F 3 (p,Q) = F 3 (K v ,X), (124) 

by means of the equation 

E?{X t K v )= m fj£' X \ (125) 

One can find E v from (|123p , substitute it into the above and integrate with respect to K v to get 

F 3 (K V ,X) = ti In (K v +Q,), (126) 

with 



fi= ^K 2 + 2hX + 4X 2 h 2 . (127) 

Now Px, the momentum conjugate to X, can be found as 

dF 3 (K v ,X) hh> 

Px = ax — --Q(K v + ny (128) 
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The above equation will be the Hamiltonian constraint after writing K v in terms of X, f and Pf . To do this, we find 
E v from (|123p , substitute it in (|122p , find K v from vanishing of it and substitute the resulting K v into (|128p . This 
way, our total Hamiltonian will be 

H„ - ffft + H tm . J» (ft + Wf) ) W ^, t m ;| J - (129) 

The next step is to introduce the second gauge condition as 

C2 = X - b{x, t) « 0, (130) 

with 6(x, t) a function of the coordinates. Since in (|129[) Px only appears in the first term in the parenthesis, the 
preservation of the above constraint 

C2 = {c 2 (x) ,j dyH tot (y)] + ^^^0 (131) 

gives 

N = b. (132) 

Because / and Pf commute with Px, the evolution equations will be 



/ = U(x),J dyH tot {y)\ = U(x),J dyH trnc (y)<> , (133) 

Pf = (Pf(x), J dyH tot (y)} = ^P f (x), J dyH tme (y)} > (134) 

with the true Hamiltonian being 

Htme = b {mPf){KM,Pf) + w,Pf))J ' (135) 

The true Hamiltonian f|135|) is a local Hamiltonian density in the sense that TV has not been given in terms of an 
integral of canonical variables. This happens thanks to the present method of gauge fixing which gives an algebraic 
equation for lapse instead of a differential equation. If N was given in terms of an integral of canonical variables, then 
the Hamiltonian would have been given in terms of an integral of an integral and would have been nonlocal in that 
sense. 

This local true Hamiltonian gives the correct equations of motion for our system as it can be checked by comparing 
the equations f|133[) and (|134[) with the equations of motion derived before gauge fixing in section (jVIII B|) . and then 
substituting in them the gauge fixing conditions. 

XI. COMPARING WITH THE ORIGINAL LAGRANGIAN THEORY 

In order to make a connection with the original Lagrangian formulation of the CGHS and also to check the 
consistency of our formulation, we transform our equations of motion into the ones in the null coordinates in the 
conformal gauge and compare them to the original Lagrangian ones. 

A. Metric and other variables in null coordinates 

The first step to make a connection between the two formalisms is finding the relations between the form of the 
metric and other canonical variables in both formulations. The coordinates used in the original CGHS formulation 
are the null coordinates 

x + =X° + x l , (136) 

x-=x°-x\ (137) 
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where x° — t and x 1 — x are the coordinates used here up to now. We can find the relation between the components 
of the null and non-null metrics using the general transformation 

dx a dx 
9ab = ^^~ 9a ' b '- (138) 

In the original CGHS model, the metric components in the conformal gauge are 

9+- = -\? P , (139) 

g— = .9++ - 0. (140) 
Thus the relations between the components in two coordinate systems are 

.goo =2.g+- = -e 2 P. (141) 

g u = - 2.9+- = e 2 P. (142) 

ffoi =9io = .9++ - .9 — = 0. (143) 
The relations between the partial derivatives in the two coordinates thus become 

d t =d++d_, (144) 

d x =d+-d_, (145) 

d t d t =d+d+ + d-d- + 2d+d-, (146) 

d x d x =d+d+ + d-d- - 2d+d-, (147) 

d t d x =d+d+-d-d-. (148) 

We can also write E x in terms of the dilaton field. Using the relation between $ and <f> (the dilaton field in original 
CGHS paper) 

$ = 2\/2e-^. (149) 

and substituting it into the jjl constraint (|52p and equating this second class constraint strongly to zero, we get 

E x = i$ 2 = 2e~ 2 ^. (150) 

The variable E v can also be written as 

\P\ 



E v =^=VQ = Vm = V9E= V=2g^=Ve 2 P = e?, (151) 



where we have used the canonical transformation Q34[) and the fact that q a b has only one independent component so 
that q — gu. The details of these calculations can be found in [18j. Using the ADM formalism, the form of the metric 
components in 2D can be expressed as 

g 00 =-N 2 + qu(N 1 ) 2 , (152) 

.9ii =9n, (153) 

g Q i=-qnN\ (154) 

From this, one can find N and N 1 in terms of metric components as 

iV 1 =--- = -^-, (155) 
9ii 5ii 



'9li „ 9oi 



N = v /9ii(^V 1 ) 2 - goo = \ ^ - .9oo = \—~ ffoo- (156) 

V 9n V 5n 

Substituting (|141[) - (|143|) in the above two equations and using (|151[) yields 

7V 1 =0, (157) 



N = ^g^=eP = E*. (158) 

Now we are ready to compare the equations of motion. 
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B. The equations of motion in null coordinates 



In order to compare our equations of motion with the original ones in the CGHS paper, we need to transform ours 
into the second order form and then bring them into the null form. Starting from the equations for the matter field 
and its conjugate, if we find P f from ([M|l. substitute it into (J85J and then use (fi"44|) - (fj"48l) . (fl"5Tjl . (fi"57j) and (fi"58f we 
get 



d+d_f = 0. 



(159) 



Finding K v from flED and substituting it in ^ and using ([THl) - (fng]) . (fT50)) , (fTCTIl . (fTST)) and (fTSgf yields 



Vi = e-Pe- 2 *(e 2 *[{d+.f) 2 + (d_/) 2 ] + 4e 2 "A 2 - 45 2 - 49 2 

- 8<9 + <9_0 + 16<9+</-<9_0 + 89 +/ oa+0 + 8d_p<9_0) = 0. 



(160) 



For the next second order equation, we find K x from (fg3"f and substitute it in (|5U)) . Then upon using (|144[) - (|148)1 . 
(Unni), (HSU), (HSU), QUE]) and the values of Pf and K^ from flBH and (gl]) respectively, we get 



V*=- ^[{d + ff + (5_/) 2 ] + 2d + + 29 2 , 
+ 4d+d-p - Ad + pd+(j) - 4<9_p9_0 = 0. 



4a,a_ 



(161) 



We can follow the same procedure and find the Hamiltonian and diffeomorphism constraints in (1571) in the null 
coordinate as 



n = e~ p e-^(e 2 ^[{d+j) 2 + {d-f) 2 ] - 4e 2 "A 2 - 45 2 4> - 45 2 4> 

+ Sd+d-cj) - 169+09-0 + 8d + pd+(j) + 8d_p<9_0) = 0, 
V = e- 24> {8d + pd+4> - 8d-pd-cf>- Ad 2 + 4> + 4d 2 _cf>) + (d+f) 2 - (d-f) 2 = 0. 



(162) 
(163) 



C. Identifying our equations of motion with those of the CGHS paper 

If we compare the above equations with the equations of motion of the original CGHS model, we can note the 
following: 

The matter field equation is identically the same in both methods and is given by (|159[) . The original energy 
momentum equations, T ++ = and T__ = in [l| are combined in the diffeomorphism constraint (J163I) as 



T++ - T__ =-V 

++ 2 



e^(Ad + pd + ^~2d 2 +( j ) ) + -(d + f) 2 



-24 



(4d-pd-^-2dU) + l(d-f) 



=0. 



(164) 



The original T_| = equation can be obtained by combining the Hamiltonian constraint equation (|162|1 and equation 

(|160p as following: 



T, 



(V! - H) 



-24> 



[2d + d-4> - Ad + 4>d-4> - X 2 e 2p ) = 0, 



(165) 



and finally the dilaton field equation of motion in CGHS paper is obtained by combining V\ and V2 i.e. equations 
(fToTJ)) and (fTBT]) as 



ePe 2 



, Vi + 2^2 = -Ad+d- 



4<9 4 



2d+d-p + X 2 e 2p = 0. 



(166) 
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XII. BOUNDARY TERMS 

It is important to take care of the boundary conditions in our theory [19| . One important reason is that energy in 
general relativity is related to the surface integral or boundary term at infinity 0. 

We require that both Hamiltonian and diffeomorphism constraints 

H c = J dxNH, (167) 

D c = [ dxN l V, (168) 

be functionally differentiable. This means that if the variations of H c and D c leads to 

/oo />oo / C tt &TJ \ /*oo 

dxNSH = dx( -r^Sq 1 + —f 5 P l + / dxd x (SS H ) , (169) 

-OO J —OO \ OQ OP J J —oo 

/oo />oo /AD AD \ r°° 

dxN x 8V = / dx( ytW + TT*' + / dxdx ( SS d) > ( 17 °) 

-oo J —oo \ OQ OP / J —oo 

then for H c and D c to be functionally differentiable, we need to add —5Sh and —SSd to them respectively. This 
means that the overall variation of surface term that should be added to the variation of the action is 

"^surface = S dt (S H + S D ) , (171) 

and clearly the surface term to be added to the action will be 

^surface = - fdt(S H + S D ) . (172) 



In the language of formulation of the CGHS that has been presented so far, i.e. using the total Hamiltonian ([761) . the 
terms that obstruct functional differentiability for the diffeomorphism constraint turn out to be 

d x (5S D ) = d x [N 1 {E^8K V + P f Sf - UJE X )] (173) 

and the corresponding terms for Hamiltonian constraint are 



d x (SS H ) = d x 



V [ II&8E* + EX " SEX + 1 fHEX + EX ' fSf 
2E x2 E^ 2 E X 2Ev 2 E x E^ 2 E X 

_ 2 \ 2 SE X - EX ' 25EV - 2 EIP ' EX ' 6EX _ K * SK <p 
E^E X E^E X E x 

3 E XI2 5E X K v P f Sf 1 PfSE x 



2 Ef 2 E x2 EvE x 2 E^ 2 E X 

+a r 



o x (J^se* 

V Ef 2 E x 



(174) 



3 More precisely it is identified as the conserved quantity associated to the invariance of the action under time translations at infinity, i.e. 
under transformation generated by a timelike killing vector field at infinity. 
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So the total variation of the surface term is 

^surface = " / dt ( dxd x (5S H ) 

dt [SS H + 6S D ]™ 



dxd x {5S H ) 



- I dtl N 1 (E^SK^ + P f Sf - U X 5E X 



N 
2 



\K% 



E" 



1 f ; 



2 E x2 e^ 2 E x 2 E<<> 2 E X 
E v'e x ' 3 E xl2 1 Pf 



E^E X 2 E^ 2 E x2 2 Ef 2 E 
( E x 'f K v P f 



2\ 2 
5E X 



\Ev 2 E x EVE X 

( NE xl 
d x ( -^r^SE* 



Sf 



E X,2 SE^ K V SK V 



E^E X 



E x 



\E i P 2 E x 
If we use the prescription at infinity for the matter field 



Sf\ 



X — ^ZOO 



= 0, 



(175) 



(176) 



then the variation of the surface term will be 



SS suliacc = - fdtlN 1 {E*5K V - U X SE X ) 



N 
2 



IK 2 



E J 



1 /' : 



2 E x2 E v 2 E x 
E v»e xI 3 E x ' 2 



2-p x 



2Ef 2 E 
1 Pf 



Ev3E x 2E^ 2 E x2 2E<? 2 E X 
E X/2 SE^ K^5K, 



2X 2 



5E X 



E^E X 



E x 



NE xl 



For the present gauge fixing, one can arrive at a surface term by substituting the gauge fixings and the values 
and shift into the general form of the variation of the total surface term (|177[) . and also putting 5h(x) = = 
(since they are just functions of the coordinates): 



S ! 



rface 



= / dt 



1 6 



Ki 



h' 2 \ 



2h V v E^ 2 ) 
Comparing the above term to the standard form of the boundary term of the CGHS [20j , 



(177) 

of lapse 
Sb{x,t) 



(178) 



S 8u rf ace = -[dt[NM\™ 



oc 
oo ' 



where M is the ADM mass, and since N — b by our gauge fixing, one gets 

2h \ E^ 2 v 



By using (|123[) . the above formula can also be written as 

M =X + 2h\ 2 
=b + 2h\ 2 

where we used the second gauge fixing in the second line. 
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XIII. CONCLUSION 

We have analyzed the CGHS model without any conformal transformation in terms of new variables similar to the 
standard Ashtekar variables of the 3+1 spherically symmetric model. This means that not only the variables have 
direct geometric meaning and there is no need to turn everything back to their directly-geometric form at the end, but 
also it might be much easier to read off the physics out of the this formulation. Then, by means of rescaling lapse and 
shift, the Hamiltonian constraint of the system was cast into a form such that it commutes with itself both in vacuum 
and coupled-to-matter cases. This makes the system suitable for analysis using loop quantum gravity techniques 
because, among other things, the system is expressed in Ashtekar-like variables and the algebra of constraints is 
now a Lie algebra. In the next step, we solved both vacuum and coupled-to-matter cases completely classically by 
introducing two gauge fixings for each case. In the case where there is matter coupling, we arrive at a local true 
Hamiltonian hence leading to local equations of motion. 

As one possible future direction, in not-gauge-fixed or partially gauge fixed case, it is desirable to polymerize the 
relevant variables that captures the semiclassical behavior of loop quantum gravity. In the totally gauge fixed case, 
one can follow different quantization schemes suitable for a local true Hamiltonian theory. These can be starting 
points for addressing the issue of singularity resolution in further works. Also the possibility of completing the Dirac 
procedure is in principle there, since the constraint algebra is now a Lie algebra. In addition, it is also worth noting 
that the algebra of Hamiltonian constraint with itself is very simple and this might be a good aid in a possible process 
of quantization. 
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